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TECHNICAL  NOTE  2687 

APPLICATION  OF  TRANSONIC  SIMILARITY 
By  Adolf  Busemann 


SUMMARY 


From  a  review  of  the  different  similarity  approaches  to  compressi¬ 
ble  potential  flow,  the  meaning  and  limitations  of  transonic  similarity 
are  traced  hack  to  their  origin.  Although  the  main  text  deals  with  the 
quasi  two-dimensional  flow,  special  suggestions  for  the  case  of  axi- 
symmetrical  bodies  are  added  in  an  appendix. 


INTRODUCTION 


As  long  as  similarity  is  used  as  a  sideline  of  theory,  its  limita¬ 
tions  may  not  be  exceeded  except  for  exploratory  purposes.  As  a  tool 
of  experimental  research,  similarity  may  more  likely  vary  between  too 
narrow  and  too  broad  applications,  or  it  may  become  fixed  to  a  standard 
form  until  its  generality  is  rediscovered  from  time  to  time.  It  is  not 
so  much  the  literature  about  similarity  as  the  comprehensive  guides  that 
seem  to  be  lacking.  There  is,  of  course,  reason  for  this  deficiency, 
since  similarity,  though  based  on  a  very  simple  principle  of  classical 
mathematics,  penetrates  progressively  so  many  different  fields  in  its 
various  applications  that  the  difficulty  of  saying  enough  about  it  is 
inferior  only  to  the  difficulty  of  adding  anything  new  to  it.  In  spite 
of  this  difficulty,  this  paper  intends  to  summarize  important  points  of 
the  background  for  transonic  similarity  in  order  to  stabilize  variation 
of  opinion  concerning  its  applications.  The  discussion  is  presented,  in 
three  parts:  The  first  deals  with  the  similarity  routine,  the  second 
discusses  the  adequate  differential  equation,  and  the  third  gives  the 
resulting  hints  for  a  proper  application. 


SYMBOIS 


A  aspect  ratio 

c(x)  or  c  center-line  deviation;  prime  with  this  symbol  indicates 
differentiation  with  respect  to  x 
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CR 

CX 
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maximum  center-line  deviation 
chord-force  coefficient 
drag  coefficient 
lift  coefficient 
pitching -moment  coefficient 
normal-force  coefficient 
pressure  coefficient 
resultant-force  coefficient 
longitudinal-force  coefficient 
Mach  number 

distance  from  axis  of  revolution 
local  static  pressure 
static  pressure  at  stagnation  point 
static  pressure  of  undisturbed  stream 

radius;  prime  with  this  symbol  indicates  differentiation 
with  respect  to  x 

maximum  radius 

velocity 

specific  volume 

space  coordinates 

angle  of  attack 

meridian  angle 

ratio  of  specific  heats 

sweep  angle 


n 


0 


4 
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3 


p  density 

t  thickness  ratio 

cp  disturbance  velocity  potential  divided  by  undisturbed  stream 

velocity  (dimension,  length) 


SIMILARITY  ROUTINE 


The  similarity  routine  in  physics  is  called  dimensional  analysis 
and  is  backed  by  the  whole  philosophy  of  conceiving  dimensions.  The 
mathematical  application  of  similarity  based  on  any  chosen  differential 
equation  is  more  flexible  but  is,  on  the  other  hand,  more  problematic, 
since  the  arbitrary  element  used  in  cutting  and  trying  a  differential 
equation  does  not  keep  the  results  within  natural  limits.  The  basic 
idea  is  always  the  same:  All  the  given  and  unknown  variables  concerned 
are  changed,  in  scale  with  no  other  restriction  than  that  of  leaving 
intact  the  underlying  physical  phenomena  or  the  selected,  differential 
equation.  Every  free  change  in  scale  is  equivalent  to  one  degree  of 
freedom  of  the  similarity  transformation  or  one  more  branch  for  the 
similarity  families.  The  number  of  degrees  of  freedom  is  therefore  the 
number  of  variables  in  the  differential  equation  less  the  number  of 
correlations  which  must  be  watched  in  order  to  leave  the  laws  of  physics 
or  the  selected  differential  equation  intact.  Once  it  is  discovered 
that  every  additional  term  of  the  differential  equation  cuts  off  one 
branch  of  the  similarity  families,  a  race  for  the  shortest  differential 
equation  suitable  to  express  a  limited  class  of  interesting  phenomena 
is  the  natural  consequence.  Any  practical  computation  has,  of  course, 
the  same  benefits  of  a  shorter  equation  but  has  the  advantage  of  allowing 
tentative  negligence  of  terms  to  be  checked  during  the  calculation  itself. 

The  compressible  potential  flow  in  three-dimensional  space,  if 
expressed  by  the  disturbance  potential  of  a  parallel  flow  and  for  a  fluid 
satisfying  the  perfect-gas  laws,  has  20  terms  in  the  general  differential 
equation  formed,  out  of  6  variables:  the  Mach  number  M  of  the  basic 
flow,  the  specific-heat  ratio  7,  the  space  coordinates  x,y,z,  and  the 
disturbance  potential  cp;  thus, 
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M2)  -  (7  +  l)l^9x  -  7  +  1  M^cp  2  -  1  ~  1  M^p  2  -  — 

2  2  y 

-  (7  -  l)*6px  -  7-^  m\2  -  2L±I  M2cpy2  -  21^-1 

-  (7  -  DM2^  -  M^cpx2  -  LJL±  M^py2  -  l±l 
(cpy  +  9x9y)  -  2M\Z  (cpz  +  9X9Z)  -  2M2q?yzcpycpz  =  0 


(1) 


Because  the  physical  dimensions  are  correct,  the  differential  equa¬ 
tion  (l)  possesses  one  degree  of  freedom  for  all  potential  flows  -  the 
geometrically  similar  flow  (same  factor  for  cp,  x,  y,  and  z  with 
unchanged  M  and  7).  Before  there  is  any  hope  of  getting  more  families 
of  similar  flows,  the  number  of  terms  has  to  be  considerably  reduced. 

The  disturbance  potential  has  been  introduced  in  order  to  reduce  the 
class  of  applications  to  small-disturbance  flows.  If  that  is  done  to 
the  extreme,  the  three  first-order  terms  in  cp  would  be  the  only  ones 
to  consider.  An  improved  accuracy  requires’  inclusion  of  all  five  second- 
order  terms  together.  The  jump  from  three  to  eight  terms  goes  beyond 
the  power  of  the  similarity  for  six  variables.  It  is,  however,  not  the 
concern  of  the  similarity  routine  to  justify  any  amputated  differential 
equation.  In  order  to  show  how  the  routine  is  carried  out,  the  basic 
differential  equation  may  be  the  so-called  transonic  differential  equa¬ 
tion  containing  all  first-order  terms  and  a  favorite  second-order  term 
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xx 


(l  -  -  (7  4  l)M2cp 


xx^x 


+  9,™  +  9ZZ  =  0 


yy 


(2) 


Justification  is  contained  in  the  next  section. 

The  number  of  participating  variables  is  still  six  and  the  number 
of  terms  in  the  transonic  differential  equation  is  four,  equivalent  to 
three  ratios  of  successive  terms.  If  all  the  ratios  are  essentially 
different,  the  original  six  free  factors  are  thus  reduced  to  only  three 
independent  factors  on  account  of  the  three  ratios  that  must  remain 
unchanged..  Combinations  of  the  variables,  which  are  independent  of  the 
scale  changes  of  the  single  variables,  are  called  "invariants"  of  the 
similarity  transformation.  A  set  of  necessary  invariants  can  be  formed 
by  the  ratios  of  successive  terms  of  the  differential  equation  (2) 
(tailing  operations  according  to  scale  values)  as  follows: 
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(l  -  M2)x 

(7  +  ljM2^ 

(7  +  l)M2jAp 


=  inv 


=  inv. 


y2 


X' 


inv. 


(3) 


The  indicated  form  depends  upon  the  arrangement  of  the  terms  in 
the  differential  equation.  The  commutative  freedom  of  the  terms  is, 
however,  no  source  of  new  invariants  since  the  same  set  of  invariants 
can  "be  combined  in  innumerable  ways  into  equivalent  invariants  which 
are  constant  as  long  as  the  first  three  are  kept  constant.  This  free¬ 
dom  of  associations  is  a  necessary  part  of  the  routine.  Actually,  it 
is  one  of  the  most  valuable  properties  for  general  applications,  though 
it  is  usually  reduced  by  convention  if  applied,  to  a  special  class  of 
problems.  Such  more  conventional  invariants  1  are  the  following  simple 
combinations  of  the  original  three: 


(l  -  M2)z2 
x2 

( 1  -  M2)3z2 

(l  +  7)2mS2 


inv. 


=  inv. 


> 


y\/l  -  M2 

X 


inv. 


(b) 


After  the  freedom  of  the  similarity  transformation,  based,  on  the 
differential  equation,  is  found  to  be  satisfactory,  the  next  step  is  to 
express  the  variables  of  the  differential  equation  by  the  given  and. 
unknown  quantities  on  the  boundaries.  It  may  happen  that  the  boundary 
conditions  of  the  original  solution  transform  under  the  remaining 


^Because  the  proper  range  of  the  transonic  term  is  close  to  M  =  1 
only,  M  in  any  power  as  a  factor  may  be  dropped  and  1  -  may 
equal  2(l  -  M)  without  any  discussion. 
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freedom  in  a  manner  that  the  resulting  new  solutions  become  impractical 
in  part  or  altogether.  The  compressible  flow  in  three  dimensions  is  an 
example.  The  zero  flow  velocity  normal  to  the  body  surface  is  not  auto¬ 
matically  conserved,  and  only  for  quasi  two-dimensional  bodies  or  for 
bodies  of  revolution  with  angle  of  attack  is  it  simple  to  find  always  a 
new  body  surface  in  the  neighborhood  of  the  old  one  without  sacrificing 
one  degree  of  freedom.  (See  appendix  A.)  In  the  quasi  two-dimensional 
case  the  body  has  the  coordinates  x,  y,  Z  with  Z(x,y)  defined  by 
the  surface  requirement 


dz(x,y) 

chc 


= 


(5) 


or 


±1  ~  x  ^  T  ^  a 
x  z 


(6) 


The  coordinate  for  small  thicknesses  Z  enters  the  thickness  ratio  t 
for  the  similar  body  and,  in  like  manner,-  the  angle  of  attack  a;  whereas 
the  aspect  ratio  A  or  the  sweep  angle  A  connects  the  actual  space 
dimensions  y  and  x.  (An  elevated  control  surface,  however,  changes 
with  Z  in  thickness  and  angle  of  attack  but  with  z  in  the  position 
above  the  main  wing. )  If  the  Mach  number  M  of  the  basic  parallel  flow  <• 

and  the  specific-heat  ratio  7  of  the  gas  are  included,  the  given  quan¬ 
tities  enter  the  ’’given  invariants”  as  follows: 

(l  -  M2)3  ^ 

(r  +  i)2m4t2 

> 

A( 1  -  M2)1/2 

or  in  any  desired  combinations,  as 

Mr  -  i)1/3tAV/3 

The  unknown  quantities  create  less  trouble  with  respect  to  their 
number  than  the  given  quantities,  since  they  have  to  be  taken  one  at  a 
time.  In  flow  problems  with  given  body  shapes  and  angles,  the  unknown 
quantities  are  almost  exclusively  pressures  or  pressure  coefficients 


(7) 


(8) 
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Cp  = 


1  2 

ipU 


cp- 


(9) 


and  integrals  of  them.  The  most  important  ones  belong  to  three  classes 
of  different  powers  in  the  thickness  ratio: 


1st  class: 


do) 


2d  class.: 


CD  ~  Cc  ~  cx  ~  CpT 


(11) 


3d  class : 

..  d2cP  ,  CP  ,  z 
da  ~  da2  ~  T  ^  x 

Their  invariants  can  be  expressed  as  follows  for  a  first-class 
repr es  ent at i ve : 

Ci(  7  +  Cpd  -  M2)1/2  Cp3C7  +  1)^ 

or  -  or  - pr -  and  so  forth 


(12) 


(13) 


and,  correspondingly,  for  the  representatives  of  the  other  classes: 


Co  T 

7T-  =  inv* 

L1 


(14) 


The  result  of  the  similarity  is  simply  that,  if  the  problem  has  a 
unique  solution,  the  unknown  invariant  stays  the  same  as  long  as  all  the 
given  invariants  are  kept  the  same.  All  members  of  one  similarity  family 
have  identical  values  of  corresponding  invariants,  given  and  unknown  ones, 
including  all  possible  combinations  of  them. 
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There  is  a  simpler  way  of  expressing  this  same  idea,  which  is  the 
representation  by  symbolic  functions,  where  the  known  and  unknown 
invariants  enter  as  parameters .  The  unknown  invariant  is  generally 
split  into  the  actual  unknown  quantity  on  the  left-hand  side  and  the 
other  factors  on  the  right-hand  side  in  front  of  the  symbolic  function; 
that  is, 


(1  -  ^)3 , 

(7  +  1)2M^t2 


a\/i 


(M  /  0;  M  /  1)  (15) 


or 


Cp  = 


-2/3 


,  1/3  2/3 

(7  +  I)  M 


- ^  ~  M  ^ — ,  A (7  +  1)1//^M2/3t1/3 

(7+I )2M  t2 


(  M-»l )  ( l6 ) 


or 


(7  +  1 )M2t 
(l  -  M2)372’ 


A 


\/l  -  M2 


(M — >0  or  7— 1)  (17) 


In  this  form  there  is  the  same  freedom  of  association  in  the  presented 
parameters.  The  first  form,  for  instance,  is  convenient  for  subsonic 
flow;  the  second  fonn  is  shaped  for  convenience  at  Mach  number  1;  and 
the  third,  form  is  convenient  for  the  case  7  =  -1. 

Some  special  cases  may  illustrate  the  meaning  of  those  changes. 
The  incompressible  result  can  be  easily  obtained  from  equation  ( 17 )  by 
letting  M  equal  zero: 


Cp  =  t  h(0,A)  =  r  H( A)  (18) 

If  the  transonic  term  is  considered  to  be  unimportant,  a  choice  of  7 
equal  to  -1  is  going  to  cancel  it.  But  before  the  elimination  is 
successful  all  former  associations  have  to  be  reversed  until  only  one 
invariant  contains  (7  +  l) ;  thus,  equation  (17)  is  again  preferable 
to  equation  (l6): 
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= 


Vl  -  m2 


h(a\^T2) 


=  tAH 


iA7T 


(19) 


After  the  transonic  term  is  canceled,  some  cases  survive  an  additional 
cancellation  of  the  term  containing  (l  -  M^),  if  this  factor  is  attached 
only  to  one  invariant,  thus 


Cp  =  tA  x  Constant 

=  a  cot  A  x  Constant 


(20) 


It  results  that  some  finite  aspect  ratio  or  sweep  angle  or  finite  wind- 
tunnel  diameter  and  a  finite  thickness  ratio  or  angle  of  attack  are 
necessary  to  express  the  unknown  quantities  with  the  given  quantities. 

If  only  the  transonic  term  is  canceled,  Prandtl's  idea  of  keeping  the 
boundary- layer  development  similar  to  an  incompressible  case  gets  the 
necessary  freedom  at  subsonic  speeds.  By  requiring  Cp  to  be  invariant 

for  boundary- layer  similarity,  the  symbolic  functions  show  that  the 


thickness  ratio  <  must  be  proportional  to 


M2 


in  similar  cases: 


if 


Cp  =  inv. 


(21) 


T 

— — —  =  mv. 
\/  1  -  M2 


and. 


Av  1  -  -  inv. 

Following  this  rule  allows  predictions  to  be  made  about  the  separation 
tendencies  since  both  potential  fJ.ow  and  boundary  layer  are  kept 
similar. 
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ADEQUATE  DIFFERENTIAL  EQUATION 


For  the  application  of  the  similarity  rules  to  experiments  the 
selection  of  the  underlying  differential  equation  has  to  he  sound.  The 
preference  of  the  lower  orders  of  cp  is  based  on  the  small-disturbance 
theory,  though  some  doubts  about  practical  thickness  ratios  and  stagna¬ 
tion  regions  are  expressed  from  time  to  time.  The  selection  of  a 
convenient  second-order  term  in  preference  to  its  four  competitors  of 
the  same  level  requires  much  more  justification.  The  background  is 
the  following:  The  first-order  terms  guarantee  an  existing  solution 
of  regular  type  as  long  as  they  are  used  at  subsonic  speeds.  At  sonic 
velocity  the  first  term  disappears  and  the  two  remaining  terms  resemble 
an  incompressible  flow  in  two  dimensions 

<Pyy  +  ^zz  =  0  f22) 


This  result  is  very  simple  and  correct,  since  the  streamtubes  at  sonic 
speed  have  only  a  second-order  change  of  cross  section  with  respect  to 
pressure  changes.  (For  analytical  solution,  see  appendix  B.)  The  main 
problem  is  to  fit  the  streamtubes  around  obstacles  regardless  of  the 
pressure  (fig.  l).  The  incompressible  medium  is  represented  by  the 
cross  sections  of  the  streamtubes.  If  any  trouble  is  encountered  by 
trying  to  fit  these  streamtubes  around  given  bodies  while  measuring  the 
cross  section  perpendicular  to  the  x-axis,  a  thoughtful  pipefitter  may 
try,  on  second  thought,  to  measure  the  true  cross  section  perpendicular 
to  the  tube  axis.  This  alternative  problem  is  identical  with  the  com¬ 
plete  15-term  sonic  differential  equation  of  a  gas  with  a  straight  com¬ 
pression  law  expressed  in  coordinates  of  pressure  against  volume,  called 
the  gas  with  7  =  -1  among  the  perfect  gases  (fig.  2).  The  actual 
difference  between  both  streamtube -fitting  problems  is  visible  but  it 
does  not  have  the  significant  result  of  saving  flow  problems  from  degen¬ 
eration.  Figures  3(a)  and  3(b)  show  the  two-dimensional  flow  in  the 
linearized  and  the  complete  7  =  -1  cases;  figures  4(  a)  and  4(b)  com¬ 
pare  the  second  case  to  the-  corresponding  axisymmetrical  flow,  the 
potential  lines  being  straight  in  figure  4(a)  and  changed  to  catenaries 
in  figure  4(b).  In  all  cases  the  flow  has  to  be  terminated  by  a  jet 
boundary  before  the  pressure  differences  grow  to  infinity,  though  the 
axisymmetrical  flow  permits  much  larger  jets  at  equal  thickness  ratios 
than  the  two-dimensional  body.  In  three-dimensional  problems  of  plan 
forms  at  angles  of  attack,  pressure  differences  are  not  quite  so  often 
degenerating  and  determine  finite  total  lifts  or  even  finite  lift 
distributions . 

The  degeneration  for  all  bodies  of  finite  length  with  finite  thick¬ 
nesses  at  sonic  speeds  indicates  clearly  that  the  general  problem  is 
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beyond  the  realm  of  the  two  remaining  first-order  terms.  In  order  to 
keep  the  solutions  from  degenerating,  the  five  second-order  terms  are 
*  the  next  ones  to  take.  Four  of  the  second-order  terms  combined  with 

nine  third-order  terms  were  included  in  the  sonic  solutions  for  the  gas 
with  7  =  -1  and  did  not  save  them  from  degeneration.  The  only  second- 
order  term  having  the  crucial  factor  7  +  1  indicating  the  curvature 
of  the  compression  law  is,  therefore,  the  most  probable  one  to  bring 
about  changes.  Nevertheless,  its  preference  requires  some  better  justi¬ 
fication  than  probability  to  succeed.  Since  the  only  task  of  the  addi¬ 
tional  term  is  to  avoid  a  known  type  of  degeneration  of  the  linearized 
equation,  the  structure  of  the  flow  field  under  consideration  is  very 
special.  Strong  gradients  occur  in  the  x-direction  and  comparatively 
smaller  ones  in  the  two  other  directions,  y  and  z.  The  quintuplet 
of  second- order  terms 


K 


-(7  + 

-(7  - 
-(7  - 

-2M^cp 

-2M2Cp 


l)M2cpxxcp. 


l)M2cp  cp 

yy  x 


l)M^zz9x 


cp 

xy  y 


cp 

XZTZ 


> 


(23) 


can  therefore  be  distinguished  by  factors  above  any  finite  limit ,  so 
that  the  one  with  all  three  differentiations  in  the  x-direction  is  far 
superior  to  all  the  others,  which  have  only  one  differentiation  in 
x-direction  and  the  remaining  two  in  other  directions.  It  can  be  shown 
that  in  a  properly  limited  thickness  range  the  predominance  of  the  first 
term  within  the  whole  quintuplet  of  second-order  terms  can  be  made  larger 
than  any  desired  factor.  In  this  situation  the  so-called  transonic 
differential  equation  is  the  product  of  theoretical  reasoning  alone. 


APPLICATION 


„  By  a  combination  of  the  information  given  about  the  similarity 

routine  and  the  justification  of  the  underlying  differential  equation, 
it  follows  that  the  transonic  term  should  only  be  used  when  the  problem 
5  forces  use  of  this  step.  This  result  means  that  one  has  to  learn  at 
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first  all  those  questions  which  can  he  answered  without  the  second-order 
term.  If  the  transonic  term  is  necessary,  the  thickness  or  the  angle 
of  attack  enters  the  results  with  peculiar  fractional  powers.  Transonic 
similarity  only  correlates  as  one  family  those  conditions  that  are 
connected  by  the  almost  vertical  structural  lines  indicated,  for  con¬ 
venience,  between  4-  and  6-percent  thickness  in  figure  5*  Similarity 
is  concerned  only  with  potential  flow  of  given  contours  composed  by  the 
solid  body  including  perhaps  a  certain  type  of  separation.  Changes  in 
separation  make  the  similarity  inapplicable,  but  there  may  be  two 
favorable  classes:  one  with  negligible  separation  and  a  second  where  the 
transonic  pressure  discontinuities  locate  the  boundary- layer  separation 
inalterably.  A  connection  between  separation  points  and  discontinuities 
has  a  known  analogy  at  sharp  corners  in  incompressible  flow.  The  tran¬ 
sonic  term  cancels  the  second  degree  of  freedom  which  enabled  Prandtl 
to  keep  the  potential  flow  and  the  boundary- layer  development  similar 
at  the  same  time  in  the  subsonic  and  incompressible  regions. 


CONCLUDING  REMARKS 


All  transonic  problems  can  be  basically  understood  by  the  addition 
of  just  one  second-order  term  in  cases  where  the  linearized  approach 
fails.  The  idea  of  the  additional  term  is  not  to  improve  the  accuracy 
of  existing  linearized  results  but  to  keep  the  linearized  results  from 
degenerating.  Adding  this  term  as  the  last  resort  in  the  similarity 
approach  shrinks  the  freedom  to  thickness  families  with  different 
boundary- layer  development.  An  experimental  check  is  not  easily  possible, 
since  below  about  4-percent  thickness  the  experimental  equipment  has  to 
be  beyond  the  customary  accuracy  and  above  about  6-percent  thickness 
the  transonic  term  is  not  sufficient.  The  original  5/3  power  of  the  drag 
coefficient  with  thickness  has  to  be  changed  sooner  or  later  to  a  first- 
power  law  because  of  the  limited  positive  and  negative  pressures  at 
sonic  speed.  There  is,  however,  a  very  good  theoretical  reason  for  taking 
transonic  similarity  as  granted  within  the  proper  limits  and  as  a  means 
of  extrapolation  with  caution  toward  larger  thicknesses.  As  soon  as 
the  linearized  theory  is  able  to  work  reasonably,  the  transonic  approach 
would  be  unwise  and  unjustified.  At  that  time  the  possibility  of 
including  the  boundary- layer  development  in  the  similarity  routine  is 
present;  the  peculiar  powers  of  the  thickness  or  angle  of  attack  dis¬ 
appear;  and  the  emphasis  can  be  placed  on  the  proper  variables. 
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APPENDIX  A 
BOUNDARY  CONDITIONS 


When  the  differential  equation  deals  with  potentials  hut  the 
boundary  conditions  are  streamlines,  it  cannot  always  be  expected  that 
similar  potential  fields  fit  to  similar  boundary  conditions.  For 
instance  the  two-dimensional  flow  may  stretch  the  potential  field  with 
Mach  number  but  shrink  the  body  thickness  perpendicular  to  the  main 
flow  direction.  The  observable  pressures  in  such  cases  on  the  surface 
of  the  bodies  shift  relative  to  the  potential  field,  but  the  pressure 
changes  at  small  disturbance  bodies  due  to  these  small  shifts  are  negli¬ 
gible  in  two-dimensional  cases,  and  the  pressures  may  be  called  the 
pressures  at  the  middle  plane  of  the  thin  wing.  In  three-dimensional 
flow  such  an  approximation  is  not  generally  permissible  as  may  be  demon¬ 
strated  for  the  slender  body  with  circular  sections. 

If  a  body  is  given  by  the  radii  of  its  circular  sections  R(x) 
and  a  center  line  c(x)  =  .z,  the  potential  distribution  in  its  neighbor¬ 
hood  may  be  arranged  according  to  the  linearized  or  parabolic  differen¬ 
tial  equation  at  sonic  speeds  of  appendix  B.  This  potential  requires 
h  just  one  more  term  than  the  proper  axis ymmetri cal  solution  (equation  (B4) ) 

because  of  the  yawed  or  curved  center  line  c(x)  and  is  written 


cp(x,y,z)  =  F0(x)  +  RR*  logN/y2  +  (z  -  c)2  -  c'r2 


z  -  c 


2  /  \  2 
y  +  (z  -  c) 


(Al) 


Gradients  in  the  y-  and  z-  directions  are 


_t-i  *  y  t„o  2y(  z  -  c) 

<Py  =  m  -5 - T - ^  +  c  R 


y  +  (z  -  c)' 


jy2  +  (z  -  c^2 


(A2) 


9-7  =  RR 


Z  -  C 


1  2  (z  -  c)2  -  y2 
+  c  R  —A - - J 


£  +  (z  -  c)^  [y2  +  ( z  -  c  )0  2 


(A3) 


J 
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These  gradients  add  to  the  radial  velocity  with  the  following  factors: 


Tr  =  9y  ■  .  7  +  TZ  ,  — 

Vy2  +  (z  -  c)2  \/y2  +  (z  -  c)2 


This  expression  shows  that  the  second  and  third  terms  of  equation  (Al) 
take  care  of  the  growth  in  cross  section  R*  and  the  slope  of  the 
center  line  c  *  as  required  for  a  solid  body  contour.  The  first  term 
does  not  enter  the  boundary  conditions  nor  does  the  arbitrary  constant 
of  the  logarithm.  Both  these  free  elements  are  connected,  since  Fq(x) 
can  be  combined  with  any  product  KR*  x  Constant.  The  y  and  z  c 
values  in  the  logarithm  may,  therefore,  be  considered  as  expressed  in 
the  proper  scale  for  the  chosen  similarity.  If  that  is  done,  the  ratios 
of  the  three  terms  should  be  invariant  for  corresponding  points.  It 
is  easy  to  see  that  the  center-line  deviation  c  has  to  follow  y 
and.  z  in  their  scales  and  that  the  first  term  in  equation  (Al)  repre¬ 
sents  one  of  the  symbolic  functions  (equations  (15)  to  ( 17 ) ) ,  if  written 
as  follows : 


F0(x)  =  f0^,  invariants^  or  ^2—  g^i,  invariants^  (A5) 

Rq  and  x0  being  the  maximum  radius  and  length  of  the  circular  body, 
respectively.  The  new  function  Tq  with  only  one  geometrical  variable 
near  the  axis  stays  the  same  as  long  as  the  similarity  is  conserved. 

It  enters  the  pressure  distribution  at  the  body  surface,  but  it  cannot 
create  lift.  The  invariant  form  of  the  second  term  depends  on  the 
chosen  differential  equation  for  the  complete  flow.  If  the  linearized 
flow  outside  the  transonic  range  is  preferred,  the  potential  reads 


NACA  TN  2687 


15 


,  ,  Rq2  /x  \  ,  ,  /Jy2  +  (z  -  c)2j  ( 1  -  M2) 

.p(X,y,Z)  =  j2-  r0(j-)  *  bb  logl/t - -jJ - 


.  2  z  -  c 

c  Rc  - 


y2  +  (z  -  c) 


(A6) 


In  the  transonic  range  it  is  more  appropriate  to  use 


cp(x,y,z)  =  -2-  gQ 
*o  u 


(5) 


+  RR 


c  rR2 


y2  +  ( z  -  c)2 


(A7) 


Another  slight  difficulty  for  circular  bodies  is  caused  by  the  difference 
in  magnitude  of  the  three  components  of  the  potential  gradient  near  the 
body.  The  axial  gradient  is  so  small  that  the  squares  of  the  crossvise 
gradients  contribute  noticeably  to  the  first  order  of  the  pressure 
changes  (in  the  same  manner  as  the  time  and  space  derivatives  enter  the 
incompressible  nonsteady  Bernoulli  equation): 


Cp  =  -2cpx  -  cpy2  -  cpz2  (A8) 

On  the  body  surface  the  simplification  may  be  used. 


y  =  R  sin  5 


) 


z  -  c  =  R  cos  5 


(A9) 


With  these  notations  the  pressure  coefficient  outside  the  transonic 
range  can  be  expressed  by  the  symbolic  function  f  and  some  extra 
terms : 


1 6 


NACA  TN  2687 


x_  cq2(1  -  M2) 
Xo2 


4. 


+ 


(2Rc"  +  J+R fc  * ) cos  5  +  2cr^cos  26  (AlO) 

In  the  transonic  range  the  corresponding  function  g  is  connected  with 
slightly  adjusted  extra  terms: 


C 


P 


x  x02(l  -  m2) 
Xo' 

Rq2M 2(7  +  1) 


co2R02M2(7  +  1) 


2(RR"  +  R'2)  log 


rrqm\/7 


+  1 


-  R 


,2 


(2RcM  +  ta'c^cos  6  +  2c*^cos  25  (All) 

It  is  obvious  that  the  extra  terms  can  only  be  used  as  long  as  the 
second  derivatives  of  R  and  c  are  limited  wherever  R  is  not  zero. 
Equations  (AlO)  and  (All),  therefore,  give  proper  information  only  if 
R  is  zero  on  both  ends  of  a  body  of  finite  length  or  the  supersonic 
character  of  the  flow  cancels  upstream  effects.  If,  however,  a  stream¬ 
lined  body,  having  limited  second-order  derivatives  throughout,  is  chosen, 
the  extra  terms  create  neither  lift  nor  drag  but  moment: 


L  =  q 


dx  dy  =  -  2jtq 


(A12) 


D  =  q 


Cp  dy  dx  =  q 


R°2  „ 

— -  g  dy  dz  +  n 

xo 


2R2R ' 2  log^i  +  R2c  ,2 


R 


x0 


_I0 

(A13) 


Mp  =  -q 


x  dx  dy  =  2jtq 


(A14) 


Since  the  first  term  creates  no  lift  either,  the  lift  and  extra  drag  are 
features  of  unstreamlined.  bodies  with  blunt  bases  or  other  disconti¬ 
nuities.  The  term  of  the  symbolic  functions  containing  the  maximum 
center  deviation  c0  may  be  of  negligible  effect  and  therefore  omitted. 
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APPENDIX  B 

SOLUTIONS  OF  THE  PARABOLIC  EQUATION 
The  sonic  differential  equation  in  x,y,z 


<Pyy  +  ?zz  =  0 


(Bl) 


is  of  the  parabolic  type;  there  is  no  direct  relation  to  changes  in 
the  main  flaw  direction  x.  For  a  real  variable  x  and  a  complex  vari¬ 
able  y  +  iz  and  its  conjugate  y  -  iz,  the  general  solution  may  be 
written 


cp(x,y,z)  =  f-j^x,  y  +  iz)  +  fg(x,  y  -  iz) 
=  R.P.  [f(x,  y  +  izTJ 


(B2) 


The  notation  R.P.  (...)  is  used  to  identify  the  real  part  of  a  complex 
number.  If  this  general  solution  is  specialized  to  two-dimensional 
flow  in  x  and  y,  only  two  powers  of  the  complex  number  y  +  iz  are 
functions  not  containing  z  in  the  real  part;  they  are  (y  +  iz)1"* 
and  (y  +  iz)l.  This  result  leads  to  the  two-dimensional  solution: 


<P(x,y) 


R.P. 


(y  +  iz )°  +  fq(x) 


(y  + 


=  f0(x)  +  fi(x)  y 


(B3). 


If  an  axisymmetrical  flow  is  wanted,  there  are  again  two  functions 
whose  real  parts  have  the  proper  combination  +  z2  =  r2#  They  are 
(y  +  iz)0  and.  log(y  +  iz).  The  corresponding  solution  is 


cp(x,r)  =  R.P.  gix)  (y  +  iz)°  +  F-^x)  log(y  +  iz^j 


=  fq(x)  +  Fl(x)  log  r 


(B4) 


¥ 
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It  is  obvious  that  neither  of  these  solutions  stays  finite  "between  the 
body  at  small  values  of  y  or  r  and  infinity.  They  are  useful  if 
applied  to  a  thin  or  slender  "body  (for  the  geometry  of  the  streamlines) 
in  a  finite  jet  with  boundaries  at  y*  or  r.  (for  the  pressure 

distribution)  as  demonstrated  in  figures  3  and  4. 


Figure  1.-  Sonic  streamtube 


Figure  2.-  Compression  law 
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Figure  3*-  Two-dimensional  sonic  flow. 


1 


* 


(a)  Two-dimensional. 


(b)  Axi symmetrical. 


♦ 


* 


Figure  4.-  Sonic  flow. 
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